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HORIZONTAL  IRRADIANCE  MODEL  FOR  SELF-CONSISTENT 
RADIATIVE  TRANSFER  APPROACH  IN  NATURAL  WATER 
WITH  ARBITRARY  VALUES  OF  INHERENT 
OPTICAL  PROPERTIES 
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1.  INTRODUCTION 

The  purpose  of  this  paper  is  to  obtain  tbe  simplest  analytic  solution  to  the  radiative  transfer  problem  in  a 
shallow  water  with  wavy  surface  and  reflecting  bottom  that  takes  into  account  all  orders  of  scattering  inside 
water  as  well  as  all  orders  of  scattering  (or  reflection)  from  bottom  and  sea  surface.  The  approach  used  here  is  a 
self-consistent  approach  (SCA)  [1-3]  that  was  originally  devised  to  obtain  upward  and  downward  irradiances  in 
shallow  sea  with  arbitrary  inherent  optical  properties.  Potentially  SCA  is  capable  to  produce  radiances  of  total 
light  distribution  in  shallow  sea  because  the  main  SCA  radiative  transfer  equation  has  a  source  term  that  depends 
only  on  radiance  of  a  source  light  and  irradiances  of  scattered  light.  This  potentiality  was  utilized  in  this  paper  to 
obtain  total  radiances  of  light  inside  shallow  sea  and  horizontal  irradiances  that  can  be  used  for  estimation  of 
horizontal  visibility.  Calculation  of  horizontal  visibility  is  considered  to  be  the  next  task  and  can  be  easily 
performed  using  the  results  obtained  in  this  paper. 

2.  RADIATIVE  TRANSFER  PROBLEM  IN  SHALLOW  WATER  WITH  COMBINED 
ILLUMINATION  BY  SKY  AND  SUN 

We  start  from  the  standard  scalar  radiative  transfer  equation  for  total  radiances  of  light  1^  in  the  ocean 

depth: 

(  d  \  b  f2;r  1*1 

+  =  d(p']_d^'p(^cosY)L,{z,n',(p'),  (l) 

here  c  =  a  +  b  is  an  attenuation  (extinction)  coefficient,  a  is  an  absorption  coefficient,  b  is  a  beam  scattering 
coefficient,  p(cosy)  is  a  scattering  phase  function  of  seawater  normalized  as  0.5  J  ^p(cos)')siny(5fy=  1, 

)'  is  a  scattering  angle,  z  is  a  Cartesian  depth  coordinate  originated  in  the  water  surface  (z  =  0)  and  going 
down  to  the  bottom  ( z  =  Zg ).  the  .x  —  _y  plane  coincides  with  the  water  surface  when  it  is  calm;  the  zenith  angle 

G=  cos  ^  H  is  measured  from  Oz  axis,  and  the  azimuth  angle  (p  is  measured  from  Ox  axis.  The  scattering 
angle  y 's  determined  by  the  following  formula: 


here  fj.=  cosd,  p'  =  cos6' .  The  source  portion  of  light  just  below  the  sea  surface  consist  of  the  following 
two  components:  1)  diffuse  component  originated  from  the  light  of  the  sky  transformed  by  the  wavy  surface  into 
Lambertian  source,  and  2)  direct  sunlight  penetrated  into  the  sea  water. 

We  consider  the  case  of  shallow  coastal  part  of  the  ocean  with  the  bottom  located  at  depth  Zg  and 
reflecting  light  Lambertially  with  bottom  albedo  Ag .  The  depth  of  the  ocean  may  be  less  than  one  optical  depth, 
1/c. 

According  to  previous  investigations  [4]  the  wavy  sea  surface  in  the  first  approximation  may  be 
considered  as  a  Lambertian  reflector.  The  albedo  of  the  surface  Ag  for  ascending  light  is  a  function  of  wind 
speed,  sea  foam  albedo,  and  refraction  coefficient  of  seawater. 

We  solve  Eq.  (1)  for  in-water  radiances  using  self-consistent  approach  [1]  to  the  radiative  transfer  of  light 
in  seawater.  This  approach  consists  of  reducing  the  scattering  phase  function  p(cos  /)  to  the  one-parameter 
transport  phase  function,  and  applying  two  experimentally-based  conditions  that  connect  the  average  cosines  of 
diffuse  portion  of  light: 

p{cosy)^  Pj(^cosy)=2B  +  2{l-2B)d(l-cosy),  d(}-cosy)=27t6{(p-(p')5iiJ.-j/),  (3) 
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=  /l„  =  l/(2  +  /f),  (4) 

fO  _  _  ' 

here  B  =  0.5  J  ^  p{pL)dlJL  is  a  backscattering  probability,  is  a  downward  average  cosine,  is  an  upward 

average  cosine,  and  Ju  is  a  total  average  cosine  [1-3].  The  average  cosines  will  be  defined  through  irradiances  of 
light  later.  Here  and  further  5(x)  denotes  a  Dirack’s  delta. 

The  self-consistent  approach  has  one  huge  advantage  over  all  other  approximate  theories:  it  is  valid  for  all 
possible  values  of  inherent  optical  properties:  extinction  c,  scattering  b,  and  absorption,  coefficients,  i.e.  it  is 
good  for  0  <  a,b,c  <  °° . 

By  applying  a  self-consistent  condition  to  Eq.  (1)  we  obtain  the  following  equation: 

(  d  \  bo  f  f  * 

\li-^  +  ajL,iz,fi,(p)=^]^  d(p'j  d^'L,(z,n',(p'),  (5) 

where  OC  =  a  +  2bg  is  a  re-normalized  extinction  coefficient,  and  bg  =  bB  is  a  backscattering  coefficient. 

Now,  let  us  split  the  total  light  radiance  L,  into  two  components:  direct  and  scattered  L  light, 
L^  —  L^  +  L.li  the  penetrated  light  of  the  sky  and  sun  just  below  of  the  sea  surface  is  expressed  as: 

(0,fi,(p)  =  Loo  +  Lgg  8(n  -  fig)  5i(p  -  n) ,  (6) 


where  Lq^  is  a  radiance  of  diffuse  light,  Lpj  is  a  radiance  of  sunlight,  then  the  radiance  of  the  source  will  be: 


LAz,fl,(p)  = 


+Lgse  “^^S(jl-flg)8((p-7t),  fl>0, 

0,  fi<0, 


here  fig  =  cos^q  is  a  cosine  of  the  angle  6g  between  direction  of  direct  solar  light  in  the  water  and  axis  Oz .  It 
is  determined  through  seawater  refractive  index  and  solar  zenith  angle  z^  via  Snellius  law, 
sinzg,  =  n„sin0o,or 


By  inserting  Li(z,fl,(p)=  L^(z,fl,(p)+  L{z,fl)  into  Eq.  (5)  with  given  by  Eq.  (7),  we  have  the  following 
equation  for  scattered  radiance  of  light  inside  seawater: 


'd\  f> 

fi^  +  a  L{z,fi)=bg\  dfi'L{z,fi')  +  Q{z), 

,  dz  y 


q{r)=]'gdfl 


is  the  source  function. 

To  proceed  further  from  Eqs.  (9)-(10)  let  us  define  diffuse  irradiances  and  average  cosines  as  follows: 
Downward  ( )  and  upward  ( )  scalar  irradiances: 


HAz)=\'^  d(p\'^dfiL{z,fi)=2niyfiL{z,fi),  (11) 

r  2;r  r  0  rO 

=  dip]  dfiL{z,fi)=2K]  dfiU,z,fi),  (12) 

Downward  ( )  and  upward  ( )  irradiances: 

d(p\'^fidfiL{j.,fi)  ^2nf^fidfiL{z,fi),  (13) 

=  “C  d(pl[^fidfiL{z,fi)=-2n\[^fidfiL{z,fi)  (14) 

fi,  =  EjH,,  K  =  EJH^,  fI  =  {E,-EAl{H,  +  H,).  (15) 

Using  Eqs.  (1 1)-(15)  we  can  rewrite  Eq.  (9)  as  follows: 
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+  (^L{z,^)^^[{2-Ji)E,{z)  +  {l  +  mM]  +  Q{z) ,  (16) 

where  radianee  of  scattered  light  is  defined  through  the  irradiances  of  that  light  and  source  function  of  external 
light  given  by  Eq.  (10). 

Equation  (16)  can  be  easily  solved  if  we  know  analytical  expressions  for  downward  and  upward 
irradiances  and  connection  between  average  cosine  ju  and  inherent  optical  properties  a  and  bg.  By  integrating 
Eq.  (16)  over  downward  and  upward  parts  of  a  solid  angle,  we  obtain  the  following  system  of  equations  for 
downward  and  upward  irradiances: 

Likiz)E^(z)=  fiiz),  1=1,2,  (17) 

where  index  1  corresponds  to  a  subscript  d,  index  2  corresponds  to  a  subscript  u,  a  summation  is  assumed  over 
repeated  indices,  fi(z)  =  'fiiz)  =  /(z)  =  27t Q(z)  ,  and  the  differential  matrix  operator  is  given  by  the 
following  equation 

d  _ 

—  +  {2-ll){a  +  bg)  -{2  +  ii)bg 

W=  ^  ^  ,  (18) 

—--  +  {2  +  ii)(a  +  bg) 

dz 

with  the  average  cosine  JI  defined  through  the  inherent  optical  properties  a  and  bg  as  follows: 


where  g  =  bg/ia  +  bg)  is  a  Gordon’s  parameter. 

In  order  to  solve  Eqs.  (17)  and,  subsequently,  Eq.  (16)  we  need  to  formulate  boundary  conditions. 
Because  we  assume  that  both  upper  (sea  surface)  and  lower  (sea  bottom)  boundaries  are  Lambertian  reflectors, 
we  define  the  conditions  as  follows: 

Eu{Zb)=  Ag[Ej{zg)+El{zB)l  £:/0)  =  ^j£„(0),  (20) 

where 

KiZg)  =  d(pl[iidiiL^{Zg,li,(p)=  2;r  +  Lo,  ,  (21) 

is  a  downward  irradiance  by  the  light  of  the  source  at  the  bottom.  The  boundary  conditions  (20)-(21)  will  be 
used  for  Eqs.  (17)  as  well  as  Eq.  (17). 

3.  SOLUTION  OF  RADIATIVE  TRANSFER  PROBLEM  FOR  VERTICAL  IRRADIANCES  OF  LIGHT 


Let  us  solve  Eqs.  (17)  for  vertical  (downward  and  upward)  irradiances  of  underwater  diffuse  light  with  the 
boundary  conditions  given  by  Eqs.  (20)-(21).  The  complete  solution  to  this  problem  is  the  sum  of  general  and 
partial  solutions  to  Eqs.  (17).  The  general  solution  is  a  solution  to  Eqs.  (17)  with  the  right  side  equal  to  zero.  The 
partial  solution  is  expressed  through  a  Green’s  matrix  of  Eqs.  (17).  The  Green’s  matrix  Gji^(z)  is  a  solution  to 
the  following  matrix  differential  equation: 

Lik(z)G^,(z)  =  5(,z)5g,  1=1,2,  (22) 

where  dg  is  a  2x2  unity  matrix  (or  Kronecker’s  symbol).  It  is  easy  to  show  [1]  that  the  Green’s  matrix  is 
expressed  as  follows: 


G,(z)  = 


^0  Hiz)e 

RoR^  (1-^0 


K 


Ro 

1  1(1- W 


•  l2-g  I  "• 


OC^  and  —(Xg  are  eigenvalues  of  matrix  operator  (18), 


(23) 

(24) 
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a^=]I{a  +  bg)-  A,  aQ=JI(a  +  bg)  +  A,  A  =  '^4a(a  +  2bg)+Ji^bl  ,  (25) 

and  H{x)  is  a  Heavyside’s  or  step  function: 

fl,  a:>0, 

(26) 

[0,  a:  <  0. 

Using  Green’s  function  (23)  we  can  write  the  solution  to  Eqs.  (17)  as  follows: 

+Dnie°‘’''  +  j‘ Gi^(z-z)Mz')dz',i=l,2,  m,  =  ,  n,=  “  .  (27) 

r~ll  ^  II 

By  inserting  Eq.  (23)  into  Eq.  (27)  and  introducing  the  following  three  auxilary  functions, 

g(z)=27r\lQ(z')e“-^'dz' ,  h(z)=2nl^^‘‘Q(z')e'^^^'dz' ,  m(T)  =  2 jV'VrfAf ,  (28) 

we  can  simplify  Eqs.  (27)  and  obtain  the  following  equations  for  downward  and  upward  diffuse  irradiances  of 
light  in  shallow  water  body  or  arbitrary  turbidity: 

E,(z)  =  A  +  DR,  +  r^^^e-“~^g{z)  +^o  (29) 

ESz)=AR^  +0  +  (30) 

The  two  coefficients  A  and  D  in  Eqs.  (29)-(30)  are  fully  determined  through  the  two  boundary  conditions  and 
can  be  written  as  follows: 

A,E’Az.)-HA,-R.mhm+g(z.)]a,e-^-’‘ 


i4c  ~  Rn  1  "I"  ■^0 

and 

E^g(z)  =  KL,^m(az)+HoL,se~“^^  (33) 

is  a  downward  irradiance  by  external  sources  (direct  and  diffuse  solar  light  and  diffuse  light  of  the  sky 
penetrated  into  the  water  body). 

4.  SOLUTION  OF  RADIATIVE  TRANSFER  PROBLEM  FOR  RADIANCES  OF  SCATTERED  LIGHT 

Equations  (9)-(10)  for  radiances  can  be  rewritten  as: 

f  d  V,  , 


^■^  +  ajL{z,n)=Qr(z), 


QTiz)  =  bg^—[i2-Ji)E,{z)  +  (2+JI)E„{z)+L,,e-^^'>‘’>]+I^^q{az)j,  (35) 

where  Eg{z)  and  Ej(z)  are  given  by  Eqs.  (29)-(33). 

With  known  right  side  the  solution  of  Eq.  (34)  for  scattered  light  radiances  can  be  obtained  as  a  sum  of 
partial  and  general  solutions: 

f  C  if  ^  0 

^z.^^)  =  l‘dz'Q,^z')G^z-z',^l)+\'  ■  (37) 


ir'  //  >  n 

where  G{z,fl)  is  a  Green’s  function  of  Eq.  (34)  that  satisfies  the  following  equation: 


/X— +  ajG(z,/z)  =  5(z), 
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and  is  expressed  as  follows: 


l/fl  ’ 


(39) 


where  H{x)  is  a  Heavyside’s  function  defined  by  Eq.  (26).  The,  insertion  of  Eq.  (39)  into  Eq.  (37)  gives  the 
following  result, 


Q  +  jj"  fir  (zO dz'l\ /f  l],  n<  0. 


(40) 


The  constants  C,  and  Cj  are  determined  by  the  boundary  conditions  (20)-(21), 


Cj  —  (Cj  +  £fl  ) ,  C2  —  Ag 


where 


4*  ^0  4*  -^OD  d-  IIq  (Lp^  /tt)  g  ^  ^ 

m^(l-AaAjp) 

fb  =  QT(z)n{az)dz,  Eg  =2jj  Qriz)n[a(zB  -  z)]dz, 
m_  =  m(azg)  =  2j^^.dn  >0,  =  m(-azg)  =  2j^iJ.diLi  >  0. 


(41) 

(42) 

(43) 

(44) 


Equations  (40)-(44)  fully  define  scattered  light  radiances.  By  adding  direct  radiance  (7)  we  obtain  the  total 
radiances  in  the  shallow  sea: 


L,(z,ll,(p)  = 


Lod  +  C^+  jlQr(z')e“^''^ dz'lll  +  “ /^o)'5(<P  "  >  0, 


r  Q  + 1^'*  Gr  dz'l\  H I 


/f<0.  (45) 


Equation  (45)  represents  a  self-consistent  solution  to  the  radiative  transfer  problem  for  radiance  of  light  in 
shallow  water  with  wavy  Lambertian  surface  and  Lambertially  reflecting  bottom.  It  takes  into  account  multiple 
scattering  inside  water  column  and  from  the  bottom  and  the  surface,  and  it  is  valid  in  the  whole  range  of 
variability  of  inherent  optical  properties  0  <  B  <0.5)  and,  consequently,  can  be 

used  for  waters  with  arbitrary  turbidity  and  absorption.  For  modeling  purposes  inherent  optical  properties  can  be 
generated,  for  example,  using  approach  proposed  in  Ref.  [5]. 

5.  HORIZONTAL  IRRADIANCES  OF  LIGHT  IN  SHALLOW  WATERS 

The  horizontal  irradiance  in  shallow  water  is  a  function  of  depth  and  azimuth  angle  ^  between  the 
direction  of  solar  rays  and  the  direction  of  light  we  measure.  It  is  given  by  the  following  equation: 

E,(z,(p)  =  ni^[L,^+qmaz)  +  C,(^{-az)  +  jl'’  ^cclz- z'])Qr{z')dz''^ 

Po cos(7r - (p),  \n-(p\<n  12, 

\n-(p\>K  12. 


(46) 


+i 


0, 


here 


(47) 


the  function  'P(t)  has  an  infinite,  narrow  and  integrable  peak  at  T  =  0,  so  it  can  be  replaced  by  a  delta 
function: 


4  a 


(48) 
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This  significantly  simplifies  the  expression  for  horizontal  irradiance: 


E,,(z,(p)  =  n<[I^o  +  q]^(az)+C2^(-az)  +  —  Qt(z)\ 


+ 


(49) 


0, 


\7t-(p\>K  12, 

with  the  coefficients  Cjand  given  by  Eqs.  (41).  It  is  necessary  to  note  that  the  third  component  in  the  right 
part  of  Eq.  (49)  is  not  equal  to  zero  only  on  a  sunny  side  of  the  viewing  direction  \jl  —  (p\<K  12. 


6.  CONCLUSION 

We  obtained  here  a  simple  self-consistent  solution,  given  by  Eq.  (45),  to  the  radiative  transfer  problem  in 
a  shallow  sea  with  wavy  surface  and  diffusely  reflecting  bottom.  This  solution  takes  into  account  multiple 
scattering  of  light  inside  water  body  as  well  as  multiple  reflections  from  bottom  and  wavy  surface.  It  is  valid  for 
waters  with  arbitrary  absorption  and  turbidity.  The  obtained  solution  was  used  to  calculate  horizontal 
azimuthally  dependent  irradiances  in  the  shallow  water  with  reflecting  bottom  and  wavy  surface.  Expression 
(49)  for  horizontal  light  irradiances  is  a  principal  equation  that  will  allow  to  solve  the  problem  of  horizontal  (or 
diver)  visibility  of  an  object  submerged  in  a  shallow  water  with  arbitrary  inherent  optical  properties. 
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